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Using the relativistic Hartree-Fock approximation, we calculate the rates of atomic ionization by 
absorption of axions of the energies up to 100 keV and for an arbitrary value of the axion mass. 
We present numerical results for atoms used in the low radioactive background searches of dark 
matter {e.g. Ar, Ge and Xe), as well as the analytical formula which fits numerical calculations for 
the absorption cross sections and can be applied for other atoms, molecules and condensed matter 
systems. Using the cross-sections for the axio-electric effect, we derive the counting rates induced 
by solar axions and set limits on the axion coupling constants. 

PACS numbers: 14.70.Pw,95.35.+d,32.80.F 



I. INTRODUCTION 

The idea of dynamical adjustment of the electroweak 
vacuum that cancels the 0-angle of QCD [1| is perhaps the 
most natural solution to the strong CP problem. This 
mechanism inevitably leads to the conclusion about the 
existence of liffht pseudoscalar particle in the spectrum, 
called axion [3]. Breaking of axial U{1) symmetry by 
QCD anomaly gives a nonperturbative mass to axions 
w^ith 
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Here m^, — {m~^ -\- rri^ )~^ is the combination of quark 
masses, {qq) is the quark condensate and fa is the axion 
coupling scale. While original models linked fa to the 
weak scale, it was soon realized that it can be in fact 
arbitrarily large [3] , limited only by cosmological and as- 
trophysical considerations (see, e.g. [J]). 

While the mass of the QCD axion is rigidly linked to its 
coupling with the topological term G'^^.G'^^ via Eq. ([T|), 
any axion model allows for additional derivative type cou- 
plings to axial currents J^d^a/ fa of quarks and photons, 
that obey the shift symmetry of the axion interactions. 
Over the years a lot of experimental activity has been de- 
voted to detecting axions using interactions of this form. 
Some methods employ finite cosmological number den- 
sity of relic axions, while others use the idea of detect- 
ing axions that are produced in the solar interior. For 
a comprehensive review of axion-related phenomenology, 
see, e.g. Refs. [l-Q- 



A dedicated search for solar axions, such as CAST [8|, 
uses conversion of keV-energy axions into x-ray photons 
in the magnetic field. Although stringent constraints on 
the axion coupling constant have been imposed by such 
searches, only recently did they become competitive with 
the broad range of astrophysical constraints. 

An alternative way of detecting solar axions was pro- 
posed in Refs. [9| . The coupling of axions to electrons can 
lead to the atomic ionization and therefore be searched 
with high radio-purity materials in the underground ex- 
periments. Recent decade has seen a proliferation of 
such experiments, that source their main scientific mo- 
tivation in searching for the nuclear recoil from scatter- 
ing of weakly- interacting massive particles (WIMPs), a 
putative component of galactic dark matter. Many of 
these experiments are also able to detect ionization cre- 
ated by solar axions down to a relevant energy scale of a 
few keV. Some constraints on solar axions were already 
imposed by the CDMS experiment [ll|. This analysis 
was also extended to the absorption of the super- weakly 
interacting massive particles (super- WIMPs) that may 
also plausibly be a dark matter candidate |10!] . In case of 
the pseudoscalar particles, the latter possibility departs, 
of course, from the mass-coupling relation suggested by 
([T]). To make the distinction clear, we shall designate 
the solar axions as "massless" or relativistic, and refer to 
the massive keV-scale axions as super- WIMP possibility. 
The constraints on super- WIMP axions were improved 
recently in Ref. [l3j . 

Up until this year, the theory of axio-electric effect was 
using very simplistic formulae relating the cross sections 



of axion absorption to the photo-electric one [l2|, [ij]. 
Earher this year, the three of us have updated these cal- 
culation for the case of the massive axions using the rela- 
tivistic Hartree-Fock calculations [la]. In this paper, we 
calculate the axio-electric effect caused by axions of ar- 
bitrary mass, including the relativistic case. Convoluted 
with the flux of the axions emitted by the solar interior, 
these results would enable searching/setting limits on the 
models of light axions that have direct couplings to elec- 
trons. Such calculations are especially timely in light of 
several dark matter experiments have reporting the ex- 
cess of events over the expected background in the keV 
region [T^ (see also Ref. [lg|, where the annual mod- 
ulation of the keV-scale energy deposition is claimed). 
These results can be generalized to constraints on the 
emission of other light particles that couple to spin, as 
e.g. in models with additional gauge bosons coupled to 
the spins of electrons jl7| . 

The main set-up of our calculation is given in the next 
section. Section 3 presents the results for the axi-electric 
cross sections. Section 4 contains calculations of the ex- 
pected signal from the solar axion absorption, and the 
Appendix provides additional details on atomic calcula- 
tions. 



II. THEORY 

The Hamiltonian for the pseudoscalar axion a interact- 
ing with electrons can be written in two equivalent ways 
[l4| (see also Appendix) 
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(2) 
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where energy scale parameter fa parameterizes the 
strength of the interaction, rrie is electron mass, a is axion 
field, ip is electron Dirac field. 

Following our previous work we present the cross sec- 
tion of the atomic ionization by absorbing an axion in a 
form which contains a dimensionless function of the axion 
energy K{ea): 
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where eg is an energy scale (in our calculations eo=l a.u. 
— 27.21 eV, but it can also be any other energy unit), c 
is speed of light, v is the axion velocity in the laboratory 
frame, qq — 0.52918 x 10~® cm is Bohr radius, Ca is axion 
energy. The function K{ea) has no unknown parameters 
and it is to be found from numerical calculations. It can 
be presented in a form 



K{ea) = - 
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where c is a state in atomic core, ric and Kc are its princi- 
pal and angular quantum numbers, k is an angular quan- 
tum number for a state in the continuum. Summation 
over L saturates very rapidly, we cut it at Lmax — 3. We 
use relativistic Hartree-Fock method to calculate electron 
wave functions in the core and in the continuum. 

The form of the single-electron matrix element depends 
on the form of the Hamiltonian for the axion-electron in- 
teraction (see Appendix for details). The first form (see 
formula (|A9IIAll|) in the Appendix) is simple. However, 
it often leads to unstable results. This is due to strong 
cancellation between the PiQj and QiPj terms in the ra- 
dial integral. The cancellation is of the order of l/(Za)^ 
which means that the formula can be reliably used only 
for heavy atoms (e.g., Xe). 

Second form of the single-electron matrix element (see 
formulas (|A16p . (|A17p and (|A18|) in the Appendix) is 
more complicated. However it is more convenient for the 
calculations since it gives stable results. In spite of some 
numerical problems, comparing calculations with two dif- 
ferent expressions is a valuable test of the calculations. 
Two forms of the Hamiltonian must give the same re- 
sults when exact electron wave functions are used. Since 
we use the Hartree-Fock wave functions we can have only 
approximate agreement between results. Therefore, com- 
paring the results is not only a test for the computer code 
but also a test for the quality of the wave functions used. 
In our experience the results agree within 10% for the 
cases when first form gives stable answers. The term 
"stable" means that variation of the axion energy leads 
to smooth change in the absorption cross section. 

Note that all formulas in the Appendix are for a closed- 
shell atom. However, this is inessential in our case. We 
consider axion energies (ca > 1 keV) for which the ef- 
fect is strongly dominated by inner closed shells while 
contribution from open valence shells is small and can 
be neglected. This means that the results can be used 
for any atom or ion with closed inner shells. They can 
also be used for molecules and condensed matter systems 
since inner atomic states depend very little on the envi- 
ronment. 

In our present calculations the axion absorption cross 
section depends on its mass. The only expression which 
depends on axion mass explicitly is the expression for the 
axion wave vector 



1 
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ka = — Vfa- {maC^f. 

In an extreme case of heavy axion (niaC^ = Ca), axion 
wave vector ka = and only term with L — contributes 
to the summation in ([5]) . This case was considered in our 



previous work [15 1 



L^c.K 



In present work we consider both these cases. First 
case gives us a test for the computer code. The results are 
the same as in our previous calculations [l^ . Second case 
gives new results for axion absorption cross section by 
atoms. According to Ref. [14] the ratio of the absorption 
probabilities for these two extreme cases is equal to 2/3 
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FIG. 1: Dimensionless factor K (see formula ([5])) in the ion- 
ization cross sections of Ar by axion. Solid line - massive 
axion {mat? ~ £a), dotted line - massless axion {nia = 0). 



FIG. 2: As on Fig. [T] but for Ge. 



in the non-relativistic limit 

o-a("ia = 0)c 



CTa{rnaC^ = ea)v 



(7) 



Below we will discuss relativistic corrections to this for- 
mula. 



III. CALCULATIONS OF THE CROSS 
SECTIONS 

Figures I1I2I and |3] show the results of the relativistic 
Hartree-Fock calculations for Ar, Ge and Xe of the di- 
mensionless function of the axion energy K{ea) which 
stands in the expression for the cross section of atom 
ionization by axion (see formula (|3])). Many body and 
relativistic effects beyond the RHF method are ignored 
and the final electron state in the continuum is calculated 
in the same potential as initial core state. The accuracy 
of this approximation is few percents due to dominating 
contribution from the inner-most core states Is, 2s and 
2p. For these states the many-body effects are small due 
to strong nuclear field. 

Solid lines on Figures [U [5] and [3] correspond to the case 
when all axion energy is due to its mass {niaC^ = to)- 
This is the same case as was considered in are previous 
work [l5|. Dotted line corresponds to the case of the 
massless axion (ttIo = 0). One can see that the ratio of 
the cross sections is indeed close to 2/3 at low energies 
(see formula ([7])). However, the ratio becomes larger at 
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FIG. 3: As on Fig. [l] but for Xe. 



high energies. For sufficiently high axion energy the ab- 
sorption cross section for massless axion becomes larger 
than those for a massive axion. This is due to the rela- 
tivistic effects. We found a formula for the ratio of the 
cross sections which fits very well the numerical calcula- 



TABLE I: Hartree-Fock energies of the core states of Na, Ar, 
Ge, I and Xe (atomic units, 1 a.u.= 27.21 eV). 



Atom 


Na 


Ar 


Ge 


I 


Xe 


Z 


11 


18 


32 


53 


54 


lSl/2 


-40.54 


-119.1 


-411.0 


-1225. 


-1277. 


2Si/2 


-2.805 


-12.41 


-53.45 


-193.0 


-202.4 


2pi/2 


-1.522 


-9.631 


-47.33 


-180.5 


-189.6 


2p3/2 


-1.514 


-9.547 


-46.14 


-169.5 


-177.7 


3Si/2 


-0.1823 


-1.286 


-7.409 


-40.52 


-43.01 


3pi/2 




-0.5953 


-5.324 


-35.34 


-37.66 


3P3/2 




-0.5878 


-5.157 


-33.21 


-35.32 


3(^3/2 






-1.616 


-24.19 


-26.02 


3^5/2 






-1.591 


-23.75 


-25.53 


4Si/2 






-0.5687 


-7.759 


-8.430 


4Pl/2 






-0.2821 


-5.868 


-6.452 


4p3/2 






-0.2730 


-5.450 


-5.982 


4d3/2 








-2.341 


-2.711 


4^5/2 








-2.274 


-2.633 


5Si/2 








-0.8762 


-1.010 


5pi/2 








-0.4341 


-0.4925 


5p3/2 








-0.3903 


-0.4398 



tions. The formula reads 



R{.Z,ta) 



Oa{ma = 0) _ 2 
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4 ^a 
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Here Z is nuclear charge, ta is axion energy, emin is the 
energy of the deepest electron state in the core for which 
ionization is possible. Note that all states in the core 
have negative energies, therefore ionization threshold cor- 
responds to the condition ta — —£min- Hartree-Fock en- 
ergies of all core states for Na, Ar, Ge, I and Xe are 
presented in Table ID The eg parameter in ([8]) is the 
energy unit. First term on the right-hand side of ([5|) cor- 
responds to the non-relativistic limit; second term is the 
relativistic correction due to core electrons; last term is 
the relativistic correction due to the kinetic energy of the 
escaping electron. 



In our previous work [15| we presented an analytical 
formula which is an accurate fit of the numerical calcu- 
lations of the absorption cross section for the massive 
axion. The formula can be used for wide range of atoms 



and axion energies. The formula reads 



K{maC^ = ea) = Kis 
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(9) 
(10) 
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where a is the fine structure constant, Z is nuclear 
charge, e^ is axion en ergy, 62 = |e2s|, 63 = |e2p|, 
vi = v/-eis/(£is + ea), V2 = 2^-e2s/(e2s + ea), ^^3 = 
2^/— e2p/(e2p + Eq)- Here eis, e2s and e2p are the Hartree- 
Fock energies of the core states. Hartree-Fock energies of 
the Is, 2s and '2,pi/2 states of many-electron atoms can 
be found using extrapolation formulas: 



ei^ _ Z2-7.49Z + 43.39 

— {Z) - , 

eo I 

— (Z) = -0.000753^3-0.028306^2 

-0.066954Z -1-2.359052, 
-0.000739^3 _ 0.027996^2 

-F0.128526Z +1.435129. 
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The functions fi{Z) and f2{Z) in (|10lllll2p are scaling 
functions: 



/i(^,e) 



(5.368 X 10-''Z - 1.17 X 10"'')e/eo 
-0.012Z+ 1.598 (16) 

j2{Z,e) = (-1.33 X 10"^Z+ 1.17 X 10"'^)e/eo 

-0.0156Z+1.15 (17) 

To find a cross section for massless axion one should take 
formula ^ and multiply it by the factor R{Z, e^) given 
by ([5]). Therefore, for the massless axion we also have the 
results which cover the same range of atoms and energies 
as in Ref. uM. 



IV. SOLAR AXION ABSORPTION SIGNAL 

To calculate the rate of the axio-electric effect caused 
by solar axions, we first address the issue of the total ax- 
ion flux. Both continuous and line-like emission is pos- 
sible. Here we take into account the emission of solar 
axions due to their couplings to nucleons, to photons and 
to electrons. The easiest case to address is the nuclear 
case, as it leads to a characteristic Ea — 14.4 keV emis- 
sion due the nuclear transition of the '"^^Fe nucleus [la ■ 



The solar axion flux was calculated in Ref. [T^ (where 
CAST results were also used to constrain it in combina- 
tion with coupling of axions to photons) . At Earth this 
flux is given by 



2.5 
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where faN is some effective coupling constant to nucleons 
that can be related to the coupling of axions to quark 
spins. The expected counting rates of argon, germanium 
and xenon experiments are given by 
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AfafaNV/' 

where the following values for the iiT-factors are used: 

ii:Ar(14.4keV) = 329; 
XGo(14.4keV) = 2746; 
ii:xc(14.4keV) = 2930. 

These rates should provide the sensitivity to (/a/ajv)"^^^ 
in the window between 10^ and 10^ GeV. Similar strength 
constraints were derived in the recent work [20], where 
a ~ 3% annual modulation of the axion signal was ex- 
ploited in conjunction with DAMA results. (Unlike the 
signal from WIMP dark matter that is expected to have 
a maximum in June, the solar axion signal is minimized 
in early July.) We leave it to the experimental collabora- 
tions to determine the exact upper limits on solar axions 
ensuing from their results. 

If the coupling to photons is not zero, Ff^i,F^^a/{4:fa-y), 
then we can calculate the counting rate, using the axion 
flux provided in Ref. Q: 



^ = 6.02 X 1030 fiGevA%2.481g-T^ 
dea. \ fat / " 



xcm ^ S ^ 



2 
J 

kcV-1 



(22) 



Counting rate for the axio-electric effect is given by the 
product of the calculated absorption cross section and 
the flux dill). For [fafa-yf'^ normalized on 10® GeV, we 
get the counting rates plotted in Figure ID Integration 
over axion energy leads to the following total counting 
rates 
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FIG. 4: Counting rate for the axio-electric effect for Ar, Ge 
and Xe as a function of axion energy. 



Comparing this to the counting rate of the CDMS 
experiment |ll| . one can see that the equivalent of 
{fafa-jY^'^ ~ 10® GeV are being probed, as the count- 
ing rates in the window from 1.5 to 4 keV reach 
0(1 kg~^day~^keV~"'^). Similar sensitivity is achieved 
in the CoGent experiment [13l |. 

Finally, the axion flux can be created by the emission 
of the axions due to the same interaction that leads to 
atomic ionization. In this case, however, the production 
cross section is down by additional factor of El/ ml [l4| . 
and the sensitivity to fa in this case does not exceed 10^ 
GeV. 



V. CONCLUSIONS 

QCD axions represent one of the most well-motivated 
extensions of the Standard Model. Their light mass and 
small couplings allow them to be produced in the Solar 
interior and escape reaching the Earth. With the prolif- 
eration of the low-background searches of dark matter, 
one should also conduct searches of solar axions. In this 
paper we have calculated the cross sections relevant for 
these searches, improving upon the simple scaling rela- 
tions that tie the axio-electric and photo-electric effects. 

Last two years has brought a significant progress in 
sensitivity to any ionizing effects in Germanium in the 
window from 1 to 10 keV [111, llJI- Currently, the Co- 
Gent experiment has very low backgrounds in the window 
from 2 to 4 keV, where the solar axion signal is expected 
to peak. With acquiring more statistics, the sensitivity 



to the solar axions in this experiment is poised to grow. 
We also remark at this point that the excess of events 
recorded by CoGent below 1 keV does not fit the ex- 
pected shape of the spectrum from solar axions. Future 
progress in searching for solar axions may also come from 
the large scale detectors with self-shielding capabilities. 
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Appendix A: Derivation of cross-section for the 
cixio-electric effect 



There are two equivalent expressions for the La- 
grangian describing coupling of pseudoscalar axions of 
mass 771a to electrons 



Ha = 2— ^a 7/;Z757/', 
J a 



1 
Ja 



K^-- (^M«) ^7''75^. 



(Al) 
(A2) 



Here ^p is the electronic wave function and fa is a coupling 
constant. The axion field a may be represented as 

with the dispersion relation 



huj — fco 



{hc\k\r 



and N being the normalization constant. 

We will treat the axio-ionization in the independent- 
particle approximation (IPA) for atom. In the IPA, the 
atomic many-body wave function is represented by a sin- 
gle Slater determinant built from single-particle atomic 
orbitals. Then, as a result of the axio-ionization, an 
atomic electron in the initial single-particle orbital will 
be ejected into a continuum state. The standard pre- 
scription for evaluating cross-sections due to Ha and H'^ 
requires computing matrix elements of transition opera- 
tors. 



Ta 



-e*'^-""i7o75, 



TL^7;e"''■''{^k^)Jorl5■ 



(A3) 

(A4) 



In the last formula 'jokfj.j'^jB = jsko — li-a. Formal equiv- 
alence of matrix elements from pseudoscalar and axial- 
vector forms of the interactions was demonstrated in Ref. 
[l4 | with the use of the single-electron Dirac equation 
with arbitrary potential. 



Since atoms are spherically-symmetric, we employ the 
standard machinery of the angular momentum alge- 
bra [22| and use the partial wave expansion for evaluating 
cross-sections. In particular, at large values of electronic 
coordinate the continuum orbital has to go over to a sum 
of an incoming spherical and plane waves [2l| . Scattering 
wave function satisfying this boundary condition may be 
decomposed in partial waves 



Vp\ 



^N, 



P/^ 



(f^L [P) Xx) ^'-^e-^'^ w,m (r) . (A5) 



Here ri^m is a spherical spinor, Xo- is a two-component 
spinor describing spin-polarization of the photoelectron, 
the relativistic angular quantum number k = {I — j){'2j + 
1) is expressed in terms of the total j and orbital I an- 
gular momenta, and (5^ is a scattering phase shift. For 
box-normalized solutions {V is the volume of the box, p 
and E are the momentum and the energy of the electron, 
and Oifs is the fine-structure constant) 



N„ 



i2nr 
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Wave function w„m (r) may be expressed in terms of 
the large (iS^) and small (T^) components satisfying the 
radial Dirac equations 






(A6) 



For bound-state orbitals, the parameterization reads 

|n,.,777,) = - (^ ^^^(^^ ^__(„^ j . (A7) 

Axio-ionization cross-sections are proportional to 
the square of transition amplitudes. Averaging 

it over all possible spin polarizations A, mag- 
netic quantum numbers m\, and 771 and integrating 
over the directions of the ejected electron momen- 
tum, we find EA«mmt«, /l(w'pAl'^a|nbKbmb)|2(if}p = 

^pZ^Kmrnt \{wKm\Ta\nbKbmh)\^ . The same result holds 
for the averaged square of T^ operator. 

The involved matrix element {'WK.m\Ta\nbHb'mb) is be- 
tween the electronic states of definite angular momenta 
and parity. For simplifying summations over magnetic 
quantum numbers, we expand the transition operators 
into irreducible tensor operators (ITO) and then apply 
the Wigner-Eckart theorem. 

We start with the simpler case of operator (|A3p . We 
employ the conventional expansion 

e^-^"- = Y. [^] '"^^^ (^^) ^Im {k) Clm (r) , 

LM 



where Clm are the normalized spherical harmonics [22 1 
and [L] — 2L + 1. Then we reexpress the transition 
operator as 



Ta = -e^'^-"-*7075 = E *"' [^] ^LM (fc) rLAiir). (A8) 



are ITOs of rank L. A matrix element evaluated between 
two atomic orbitals reads 



LM 



Here the operators 

TLM(r) = - 17075 JL (kr) Clm (f) 



{TLM)ij = -{nimi\CLM\ - Kjirij) i / JL (kr) [Pn,Ki{r)QnjKj{r) + QniK,,{r) F„^.K^.(r)] drj 

The selection rules for matrix elements of the C-tensor require that \ji — jj\ < L < ji + jj and L + k + Ij = odd. For 
example, for L = the multipolar operator is pseudoscalar: the Too operator drives si/2 — > P1/2 transitions. Reduced 
matrix element 



(njKj||ri||njKj) = -(Ki||CL|| - Kj) ( / Jl (kr) [P„,Ki (r)(5„^.„^ (r) + (9„,„.(r) Pn^^jir)] dr\ 



(A9) 



To evaluate the cross-section we fix the coordinate system in such a way that the axion propagates along the z-axis. 
Then in Eq.(|Ml), 



Clm 1 k 



Smq and Ta = ^i^[L] tlm=o{y) 



(AlO) 



Further 



Finally, 



Yl Y. f\{wpx\Ta\nbKbmb)\^dnp=N^ ^ (2i + 1) ((K||TL||nb«:b)) 



Akto n^KfjTrif, 



n^Kbi^L 



C AtT 1 1 ■r-^ 2 

a= 2~T2~ ^ /. (2L+1) ((k TL nfeKb)) 



(All) 



n^K^jKh 



Derivation of the axio-ionization cross-section for the alternative form of the coupling H'^ is more complicated. We 
start from the multipole expansion of the T'^ operator. 



T'a = \^ {l,k, - k . a) e'*--- = ^ z^ [L] CU [k] Au^^ 



(A12) 



LM 



Because the angular dependence of this expansion is the same as in Eq. ljASp . the expression for the cross-section 
remains the same as in the Ha case, Eg. ljAlip . with the substitution tl ^ t'j^- 
The multipolar tensors t^^;/ may be derived by inverting the expansion (IA12P 



o-i 



'^z,A/(i") == -^ I dn^CLM 



{k)l^H 



sfco - k • a) e' 



Y, «'^'"'' [L'] JL' {kr) Cl'M' (r) ^ / dn^CLM (k) \i {l^k^ - k • d) Q,^' {k) 

L'M' ^ •' 



The two contributions to the integral are 



and 



— / d^kCLM [kj Q'Af [kj -i (75fco) = j-r-i (75^0) '5ll"5mm' 



— / dVLkCLM [k] CI, HI, (k) -i (-k ■ a) = --\k\i{a ■ sll'M'-lm) , 



(A13) 



(A14) 



where the components of a vector object sll'M'-.lm 

1 



{^L'M';Lm)x 

The resulting expression reads 



47r 



dn^ClM' (fc) kCLM (fc) 



V[iW\ 



■.{L'M'\Cix\LM). 



(A15) 



rLiir) = <M(r) + rs'^Mir) - Jl (fcr) Clm' (r) J* (75^0) " ^ |k| J] *^'-^+i [i'] Jl' (fcr) C^^m' (r) (a • ai,M';LA/) 



L'A/' 



where we spht the operator into the time- and space-hke contributions. Below we tabulate reduced matrix elements 

{niKi\\T'j;^\\njKj) = {niKi\\Tt'];^\\njKj) + {niKi\\Ts'j;^\\njKj), (A16) 



of the t'^j^i ITO 



{niKi\\Tt'^\\njKj) = -^koiKiWCLW - Hj) ( / Jl {kr) [Pn,Ki{r)Qn,Kj{r) - Qn,K,{r) Pn^Kjir)] drj , 



(A17) 



(UiKi TSr. n 



ir^lW ' ^L I I "'j'^3 



K,)=--M 



E 



•L'-L+l, 



k,\\Al'\\hj)[L']x 



L'=L-1.L + 1:L'>0 



f JL' (kr) [P„.,, {r)Pn,n, (r) + Q„.«, {r)Qn,^, (r) ] drj , 



with 



W\\Al'\\k) = {k'WY^Cl'M' (f) {<J ■ aL 

M' 
7/1 1/2 



M':LM 



)\\'^) 



[L,Lf/2 ^ ^ ^ ^ lo 



L' 1 L 




L' L 

1/2/' / 
1/2/ J 



(A18) 



Here the notation [Ji, ..., J„] = (2Ji + 1) ... (2 J„ + 1). The two-row quantities are the 3j-symbols and the 3x3 matrix 
in the curly brackets is the 9j-symbol. Notice that the phase i^ -^+1 entering Ea. (|A18l) is either -1-1 or -1, i.e., the 
entire expression is real. Selection rules for both time- and space-like contributions are the same as in the case of the 
tlm multipoles: \ji — jj\ < L < ji + jj and L + k + Ij = odd. 

To summarize, the cross-section for axio-ionization is given by Eq. (|Alip . with reduced matrix elements given by 
Eqs. (jA9p and (jA16[) . Notice that the derived expressions remain valid for arbitrary large values of parameter kr, i.e., 
even when the usual dipole approximation (e**''' ~ 1) breaks down. 
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